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The search for nonvolatile memory concepts has a massive impact 
on the development of nanoscopic systems with adjustable electrical 
properties. Capacitor-like structures composed of insulating mate- 
rials sandwiched between metallic electrodes are envisioned to over- 
come the limitations associated with conventional charge storage de- 
vices and may open the road to neuromorphic computing. Together 
with phase-change memories [U |2] two concepts attracted extensive 
interest [3j: resistive and magnetoresistive random access memories 
(RAM). In resistive RAM, electrochemical processes within the in- 
sulating layer enable the access to multiple resistive states |4j[5l[6l[7]- 
For magnetic RAM, ferromagnetic electrodes are used to switch the 
resistance by changing the alignment of the electrodes from parallel 
to anti-parallel and vice versa [SJ |9j I10| . Here we demonstrate that 
both effects can be observed simultaneously in nano-scale magnetic 
tunnel junctions (MTJs). The framework of a 2nd order memristive 
system is utilized to treat the tuning of electrical resistance by 
concurrent modulation of two material parameters. 



The development of nanoelectronic nonvolatile memories and neuromorphic 
computer architectures requires devices with stable resistance states and ana- 
logue capabilities to mimic, e.g., neurobiological architectures. Both concepts 
have large significance for development and improvement of functionalized ma- 
terials and nano-scaled systems with adjustable electric properties. Specifically 
engineered memristors provide controllable resistance where the definition of 
the memristor is based solely on fundamental circuit variables. The Memris- 
tor theory was formulated by Leon Chua in 1971 [12]. Chua extrapolated the 
conceptual symmetry between the resistor, inductor, and capacitor, and in- 
ferred that the memristor is a similarly fundamental device. Other scientists 
had already used fixed nonlinear flux-charge relationships, but Chua's theory 
introduces generality. The interest in this concept renewed in 2008 due to ex- 
perimental progress by Strukov et al. |13j where a memristor was realized as 
a Pt/TiO/Pt sandwich. Moreover, Strukov et al. created a link between the 
pioneering work of Chua and the field of resistive switching. The memristor "re- 
members" the history of the applied voltage v and current i since its resistivity 
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Figure 1: Simultaneous occurrence of resistive and magnetoresistive 
switching, a, A slight splitting of the i-v curve can be observed for both mag- 
netic states, which demonstrates the presence of resistive switching. Since the 
curves are highly symmetric with respect to the origin, only the first quadrant 
is shown, b, The magnetoresistive switching of the device characterized by a 
magnetic minor loop, (see also Supplementary Information Fig. 1 and 2) 



depends on the time integrals of v and i. A voltage controlled memristor can 
be described by 

u(t)=i2 ($(*))•<(*) (1) 

where $(t) = u(t')dt' is referred to as flux and R ist the (memory) resistance 
[14U15j . The memristor was generalised to memristive systems in a 1976 publi- 
cation by Chua and Kang [11] . Whereas a memristor has mathematically scalar 
state, a memristive system has vector state. The number of state variables 
is independent of, and usually greater than, the number of terminals. In this 
paper, Chua applied this model to empirically observed phenomena, including 
the Hodgkin- Huxley model of the axon and a thermistor at constant ambient 
temperature. He also described memristive systems in terms of energy storage 
and easily observed electrical characteristics. These characteristics match resis- 
tive random-access memory and phase-change memory, relating the theory to 
active areas of research. In the more general concept of an nth order memristive 
system, equation]!] turns to 

u(t) = R (w, u, t) ■ i(t) (2) 
w = f(w,u,t) (3) 

where to is a set of n state variables describing the device [11] [HI [15] . Neither 
R nor / are explicit functions of u in the case of a "pure" memristor. 

After substantial effort was invested in studies of resistive and magnetore- 
sistive systems, the next step is a combination of both effects in one device 
[TBI H"7l [18] . These systems can be treated as 2nd order memristive systems 
with w\ the magnetic state of the electrodes and W2 describing the state of 
the barrier. In this framework, we present studies on MgO based low resistive 
magnetic tunnel junctions with Co66Fe22Bi2 electrodes. The samples are char- 
acterized by tunnel magneto resistance (TMR) ratios of about 100 % and exhibit 
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Figure 2: Resistive switching (RS) of a magnetic tunnel junction. The 

RS is induced by the sinusoidal voltage plotted at the right axis. The resistance 
of both magnetic states (AP and P) is switched periodically between two resis- 
tive states (i?H and Rl). Labels are attached to the resistance curves (at V = Q) 
to clarify the definition of Ru and R^. A magnetic field of ±400 Oe provides 
a stable parallel or anti-parallel orientation of the magnetic layers which sup- 
presses spin transfer torque effects during the RS measurements. The read-out 
delays are neglected in the time scale. 

an additional bipolar resistive switching (RS) of up to 6%. The TMR ratio is 
not reduced by the operation as resistive memory and the area resistance is low 
(6.8 ri^m 2 ). Therefore, the presented results arc a leap towards a combined 
magnetoresistive and resistive memory that could be switched by spin transfer 
torque [19]. Figure 1 shows the current- voltage (i-v) curves for the parallel (P) 
and anti-parallel (AP) orientation of the Co-Fe-B electrodes. The current value 
measured for increasing voltage is slightly lower then the current measured for 
decreasing voltage. This hysteresis of the i-v curve is observed in both magnetic 
sates and is the primary characteristic of memristive behaviour. 

To further investigate this effect, we reduce the applied voltage once a second 
and measure the resistance at 20 mV with a read-out delay of 200 ms. This pro- 
cedure removes non-linear contributions resulting from the voltage-dependent 
tunnel resistance |20j and is typical for nonvolatile memories where the infor- 
mation is written at high and read at low bias. Figure 2 shows the variation 
of the MTJ's resistance induced by the application of sinusoidal voltage. The 
resistance of the AP state is switched between a high resistance (Rn) and a low 
resistance (Rl) state. A similar resistive switching is observed for the P state. 
The high and low resistance states are defined by the resistance measured at 
zero voltage, as indicated by the circles in Figure 2 and 3. Table [T] summarizes 
the results obtained from the data presented in Figure 2. 

Figure 3a shows the resistance of the P state now as a function of the applied 
voltage. Starting at i?H the resistance remains almost constant until a critical 
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i?H Rl Ai?Rg RS 



R AP 375.9 362.5 13.40 3.7% 

R P 189.6 178.9 10.70 6.0% 

AiJxMR 186.3 183.6 

TMR 98.3% 102.6% 

Table 1: Resistive (RS) and magnetoresistive (TMR) effects calculated from the 
data in Fig. 2. The index to R is denoting the (magneto)resistive state of the 
device. AR is the total resistance change due to the effect denoted by the index. 
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Figure 3: Two complementary methods to display the resistance of 
the P state. 2.5 resistance loops are plotted on top of each other, a, The 
bipolar resistive switching between two stable resistance states becomes clear 
when plotted vs. voltage, b, The hysteretic nature of resistive switching is is 
shown with clarity when plotted vs. flux. 



voltage of about 0.37 V is reached. Then, the resistance drops reaching the 
lowest value for maximum bias voltage. A slight increase of the resistance can 
be observed when the applied voltage is released. A similar switching process 
appears in opposite direction after the polarity of the applied voltage is changed 
and, again, when a critical voltage for the switching is reached (— 0.42 V). Fig- 
ure 3b shows the hysteretic nature of the resistance as a function of flux. The 
flux "applied" to the sample increases when a positive voltage is applied and de- 
creases for negative polarity. For increasing flux the resistance decreases slightly 
till a flux of about 14 Vs is reached and the resistive switching begins. The 
minimum resistance is reached at about 42 Vs before the resistance increases 
slowly reaching R^ at 84 Vs. The back-switching begins at 67 Vs leading to a 
pronounced resistive hysteresis. The resistive loops lie on top of each other, 
indicating the reproducible nature of the switching process. 

As pointed out previously [2T], a possible explanation for the RS of MTJs 
is a displacement of oxygen vacancies located at the bottom Co-Fe-B/MgO in- 
terface. We perform measurements on MTJs with an MgO barrier prepared by 
post-oxidation of an Mg layer to study this model in more detail. The bottom 
Co-Fe-B electrode is not exposed to oxygen for this preparation method and no 
oxygen vacancies at the bottom interface are expected [52]. The correspond- 
ing measurement is shown in Figure 4a. Two features are noted. First, no 
i?($)-hysteresis can be observed indicating the absence of resistive switching, 
as expected. Second, a V-shaped curve is observed with the resistance lowered 
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Figure 4: Thermal contribution to the R(&) curve, a, Reference mea- 
surement on an MTJ showing no RS. The amplitude of the applied voltage was 
350 mV, 2.5 periods are plotted on top of each other, b, The fitting-curves 
calculated by Eq.|4]are shown by dotted lines while the experimental data, the 
same as in Fig. 3b, is given by the curve in the background. The deviation 
of the experimental data from the fitting-curve due to thermal contributions is 
depicted by red lines. 




by 0.62 % in the middle of the flux-axis. At this flux the applied voltage reaches 
its maximum for both polarities. Since the tunnelling elements (and the top 
and bottom conduction lines) are heated due to the high current density the 
local temperature may not relax during the 200 ms of read-out delay. Therefore, 
a possible explanation for the reduced resistance is an elevated temperature of 
the MTJ. A similar influence of remnant heat is visible for the sample showing 
a distinct RS behaviour as denoted by the red lines in Figure 4b. 

Now, we stress the question how the presented system differs from a pure 
memristor. Despite the superposition with thermally induced resistance changes, 
the hysteresis of R(Q) means that a pure memristor can be defined at best 
"branch- wise" : 



with i? + (<i>) the resistance on the positive branch (V > 0), $g = 19.5 Vs the in- 
flection point and S + = 4.3 Vs indicating the sharpness of the switching process. 
To fit the negative branch i?~(<I>) we shift the inflection point to <I>g = 52.8 Vs 
and use 5~ = 5.3 Vs. Both fitting curves are shown in Figure 4b. However, a 
simple hysteresis as presented in Figure 3 will only result from a sinusoidal volt- 
age of constant amplitude. If the positive amplitude is alternated between two 
values, as illustrated in Figure 5, an additional i?L state is generated. In gen- 
eral, multiple resistance states both high and low can be created by according 
amplitude variation. 

In summary, simultaneous magnetoresistive and bipolar resistive switching 
of MTJs is presented. The entire structure can be understood in the framework 
of an 2nd order memristive system. Five resistance states are demonstrated 
for one MTJ proving that the combination of resistive and magnetoresistive 
switching of MTJs provides a method for multi-bit data-storage and logic. 
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Figure 5: Influence of an alternating voltage amplitude on the i?L state. 

The positive amplitude is periodically switched between 450 mV and 550 mV. 
This leads to a splitting of the low resistance state into i?Li and i?L2 separated 
by 6.7 Q. 



Methods 

The TMR systems are sputter deposited in a Singulus timaris II tool. The 
film sequence is (bottom) Ta 3/Cu-N 90/Ta 5/Pt 37 Mn 63 20/Co 70 Fe 30 2/Ru 
0.75/Co 66 Fc2 2 Bi2 2/MgO 1.3/Co 66 Fe 22 Bi 2 3/Ta 10/Cu 30/Ru 7 (top). The 
thickness is given in nm and the composition in at. %. Elliptical TMR elements 
with sizes in the range of 0.038±0.005 /im 2 are prepared by electron beam lithog- 
raphy in combination with ion beam etching (see Supplementary Information 
Fig. 3). The film sequence for the reference structure showing no RS behaviour is 
(bottom) Ta 5/Cu-N 90/Ta 5/Pt 37 Mn 63 20/Co 70 Fe 30 2.2/Ru 0.8/Co 60 Fe 2 oB 2 o 
2/MgO 0.8/native ox. ltorr 1200 s/Mg 0.3/Co 60 Fe 20 B 2 o 1.5/Ta 10/Cu 30/Ru 
7 (top). 

All transport measurements are performed using a constant voltage source at 
room temperature with the bottom electrode grounded. To induce RS we apply 
voltage pulses of up to w max = 0.65 V for one second to stress the samples by a 
current density of typically 8 x 10 6 A/cm 2 . All resistance values are measured 
at a read-out voltage of 20 mV with a read-out delay of 200 ms after each stress 
pulse. In total, the applied voltage sequence is v n (t) = v max sin(t) for every odd 
pulse-number n and v n = 20 mV for even pulse-numbers (see Supplementary 
Information Fig. 4) . 
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Figure 6: A slight splitting of the i-v curve can be observed for both magnetic 
states which demonstrates the presence of resistive switching. Since the curves 
are highly symmetric with respect to the origin, only the first quadrant is shown 
in the main article. 
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Figure 7: Resistive switching without barrier degradation. To get the presented 
data the MTJ is switched 5 times between the high resistive state i?H and the 
low resistive state i?L- The TMR ratio is obtained from an magnetic minor loop 
measured in the i?L state (open circles). Again, the MTJ is switched 5 times 
and the TMR ratio is measured, now in the i?H sate (solid circles). Then the 
voltage amplitude is increased by 250 mV and the procedure repeated. Since 
the TMR ratio remains high, a degradation of the MgO barrier due to resistive 
switching can be excluded. 
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Figure 8: Scanning electron microscope (SEM) micrographs of typical elements 
prior to etching. The width of the resist mask is 119nm x 337 nm (top) and 
130 nm x 343 nm (bottom). 
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Figure 9: Details of the measurement procedure. The applied voltage sequence 
is v n (t) = i; max sin(i) for odd pulse-numbers n and v n = 20 mV for even pulse- 
numbers. The length of the odd pulses is 1000 ms, the length of the even pulses 
is 200 ms. All resistance values are obtained after an odd (delay) pulse at 20 mV, 
whereas the resistive switching is induced by the even voltage pulses. 



12 



